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Abstract 



ff^ ' In a fundamental paper [Phys. Rev. Lett. 78, 325 (1997)] Grover showed how a quantum 

'nJ" ■ computer can find a single marked object in a database of size N by using only 0{\/N) queries of 

cn ■ 

I the oracle that identifies the object. His result was generalized to the case of finding one object in a 



subset of marked elements. We consider the following computational problem: A subset of marked 
elements is given whose number of elements is either M or K, our task is to determine which is 



P-t. the case. We show how to solve this problem with a high probability of success using iterations of 

> ' 

Ch ' Grover's basic step only, and no other algorithm. Let m be the required number of iterations; we 

! 

prove that under certain restrictions on the sizes of M and K the estimation m < Jy^^ obtains. 
This bound reproduces previous results based on more elaorate algorithms, and is known to be 
optimal up to a constant factor. Our method involves simultaneous Diophantine approximations, 
so that Grover's algorithm is conceptualized as an orbit of an ergodic automorphism of the torus. 
We comment on situations where the algorithm may be slow, and note the similarity between these 
cases and the problem of small divisors in classical mechanics. 

PACS numbers: 03.67.L.X. 
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Consider a database of = 2" elements, which are represented as the basis vectors of a 
quantum register \aj) = |a{) (8> ... (8> |a{), al E {0, 1}. The state of the register can be any 
superposition the basis vectors. Our task is to find in the database one specific element \aj). 
At our disposal is an oracle that if given the required element \aj), will mark it by rotating 
its phase by tt. Should the oracle receive a superposition of the basis elements, it will rotate 
only the branch of \aj). Grover [l| demonstrated that by using 0{y/N) calls to the oracle 
one can find the marked element \aj) with a very high probability. It was also shown P] 
that if one is asked to find any one oi K {1 < K < N) marked elements, it is possible to 



reach a high probability of success by calling the oracle 0{^jN/ K) times. 

In this paper we consider a variant of the algorithm which can solve fast the following 
problem: We know that there is a subset S of marked elements in the database, and we have 
an oracle to demarcate them. However, we do not know exactly how many elements there 
are in S", only that the number is either \S\ = M oi \S\ = K for some Q < M < K < N/2. 
We shall show that under certain restrictions on the values of M and K and their relations to 
N we can solve the problem by calling the oracle m times, where m < ^^^^]^r^ (theorem 1). 

Our result reproduces an earlier work of Nayak and Wu who obtain a solution to the 
problem with a probability of success > | after calling the oracle m < O ( -i / + ^^^m^'' ) 



times. The authors apply the counting algorithm of Brassard and his collaborators [45 , which 
involves the elaborate discrete Fourier transform |[5| on top of Grover's simpler procedure. 
By contrast, we just iterate Grover's rotation, and apply simultaneous Diophantine approxi- 
mations ^,6] to calculate the number of iterations that separate the two cases. Mathematically 
speaking, this means that we conceptualize Grover's algorithm as an orbit of a discrete dy- 
namical process on the torus T^. Note also that alHhese upper bounds are optimal, up to 



a constant factor. This was proved in |3|, see also 

First we shall briefly repeat the algorithm 
{1, 2, A^}, such that l^l = K. Let us denote 



J. 



2] of finding an element of a set S" C 
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Now, write the initial state of the register \ip) = X]i=o 1*^*) a sum of the two vectors 



in Eq O 1^) = \a)^ + J§ |/?)^, or 



1^) = cos— |a)^ + sm— 1/3)^, ^ " V 7? (2) 



Each step in Grover's algorithm transforms the present state of the register to a new- 
state G l'^'), where G is a rotation of the plane spanned by \oi)j^ and by the angle Ok- 
To perform the rotation we first call the oracle to reflect IV'') around |a)^ by introducing a 
minus sign to the \l3)j^ component; subsequently we reflect the result about After m 
iterations the state is 

1^) = cosimOK + y) + sinMK + ^) \(3) (3) 

All that is left to do is choose m that will bring as close as possible to in other 

words, wc look for an integer m which will satisfy sin(m + |)^/^ ^ 1. In case N ^ K it 
follows from the definition of 9k that m is of the order of magnitude of \J^- 

In the present problem we are given integers M < K < N, and we are told in advance 
that there is a subset S C {0, ...,A'"} of marked elements that contains either M or K 
elements but we do not know which is the case. We wish to find out whether l^"! = M or 
\S\ — K. To flnd the answer we simply apply Grover's rotations, without any change. The 
only difference from Grover's original procedure is the stopping rule, that is, the number of 
iterations required before the measurement is performed. So essentially the same algorithm 
is solving a different problem. 

We can always represent the initial state 1-0) as 

\^) = cos ^ \a)^ + sin ^ if \S\ = M, (4) 

9 9 
1^) = cos^ + sin^ 1/3)^ if \S\ = K 

With sin ^ = -y/^ and sin ^ = Our purpose is to compute the number of 

iterations m with the following property: If jS*! = M the rotation G operates in the plane 
spanned by \a)j^ and and G"^{\ip)) is close to \a)j^, the vector of — M unmarked 

elements; however, if IS"! = K the rotation G operates in the plane spanned by |q;)^ and 
\P)k while is close to the vector of K marked elements. We do not know 

in advance which is the case, but if such an integer m is found, and G has been iterated 
m times, all that is left to do is measure the quantum register. If the result is one of the 
elements of S (which we check by another query of the oracle) then it is clear with high 
probability that l^l = A' , otherwise, l^l = M. As usual, the probability of success can be 
further increased by repeating the process. 
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The existence of such an integer m follows from the theorem of Kronecker on simultaneous 
Diophantine approximations p] : let ^1,^2, ■■■,^r be irrational numbers which are linearly 
independent over the rationals, and let rii,ri2, ..■^T]r be any real numbers, and e > real. 
Then there are integers pi,p2, ...,Pr and an integer I such that 



- Vj - Pj\ < £ i = l,2,...,r (5) 
In our case r = 2, and we wish to find an odd integer I = 2m + 1 which approximates 



^1 = ^ to r/i = 0, and ^2 = |^ to 772 = | |Si- Only in rare cases such ^1 or ,^2 are rationals, 
or dependent over the rationals 9] . The trouble is that it is very hard to obtain a universal 
bound on the minimal number / that satisfy Eq In the general case of arbitrary ^j's and 
?7j's no such universal bound exists. In our more specific case, when we consider all M, K, 
and A^, it is an open problem. 

Luckily, there is an interesting range of values of M and K for which a small odd I does 
exist. Denote 7 = = ——r-^'^^-, our main result is 

' sin-i(v^)' 

Theorem a. If M < K < ^ satisfy ^/K < 16(7 — 1)'^^/N then there are natural 
numbers I and p, such that I is odd, and I < J/^nrr , and 



< 2(7-1) 



< (7 - 1) (6) 



b. Let 6 > and consider the cases where K < (1 + ^^)^M. Then the inequality (7 — 1) < | 
is satisfied, so that —p—W < ^, and also — p\ <s. 

The proof of the theorem is given towards the end of the paper. The theorem allows us 
to solve the problem with a high probability of success. Suppose that we have iterated the 
algorithm m = ^ times, m < and subsequently measured the register. Then the 

y i\ — V A'l 

probability of getting the wrong result is determined in the following way: First, suppose 
that there are K elements in S*, then the probability of getting an unmarked element after 
the measurement is 

JO 

\{G'^^\a)^\^ = cos\meK + = cos\-^) < sin2(27re) (7) 

where the last inequality follows from the theorem. Likewise, if there are M elements in IS"!, 
the probability of measuring a marked element after m iterations is 

KG-^mMf = sin^M., + ^) = sin2(^) < sm'{2ne) (8) 



again, the last inequality follows from the theorem. Note that e need not be excessively 
small. Even if we take sin^(27r£) = |, that is £ = ^, then a few repetitions of the algorithm 
will give the correct answer with overwhelming probability. 

In light of the theorem one can see Grover's algorithm as a discrete dynamical process on 
the torus T^: Consider the subset of given by D = {(/|^ (mod 1), (mod 1)) ; lodd}. If 
^ and ^ are independent over the rationals, then D is dense in T^; this is just Kronecker's 
theorem (with the slight variation that we consider only odd Ts). If Z = 1, 3, .. is taken as a 
discrete time parameter, then is a dense orbit of an ergodic dynamical system, and the 
question is how quickly it will enter a small prescribed neighborhood of (|, 0). This question 
can be generalized to more extensive searches on T'', for r > 3 (more on this below); or to 
questions concerning approximations to other points on the torus, which may be related to 
the solutions of Diophantine equations; or finally, to questions regarding continuous rather 
than discrete processes, such as adiabatic computations. We shall come back to this point 
later. 

Here are a few applications of the theorem: 

Example 1 For K — M + 1 we need m = ^ < (M + 1)A^ iterations of Grover's 
algorithm to solve the problem up to a probability of error sin^(27r£). To estimate the range 
for which this is possible note that since > ^ for < X2 < xi < 1, we have in this 

case 7 — 1 > + jj ~ 1 > Therefore, if we choose \J^^ < ^^^^ condition 

of the theorem: \fK < 16(7 — l)^\/iV is fulfilled. This means that the range of application 
of the algorithm for this case is at least M ^ -v^; and the number of steps is m 

Example 2 For K = 2M we can increase the database by adding rN artificial elements, 
out of which rM are marked, so they respond positively to the oracle. As a result we have to 
separate now between M' — (r + l)M and K' — {r + 2)M, while the total size of the database 
increases to N' — {r + 1)N. We proceed as follows: 

a. The condition that K' < (1 + -^^^M' is satisfied ifr + 1 > \ple~^ . Let r be the minimal 
integer that satisfies this inequality. 

b. Consider the new angles: 6m' = sin~"^(-^ ^^^^^^'^ ) = 9m, andOx' = ^i^^ "^(-y/^^f^)- Then 



^'-^ = €-^^ V - 1 > ^- (^^ ™^ ^^^^ "^^^^ IHS) ^ t 

for < X2 < Xi < 1; and the minimality of r from a). 

c. Consequently, if we assume < (y)^ then the condition \J~^ — ^ +1) jv < 9^^ ^ 



16(7' - 1)^ is fulfilled. 

d. Now, apply the theorem to separate between M and 2M with an error e in a number of 
steps m less than ^j^^^^f^^^-^ < 5(r + 1)^^ ~ ^^^^^ \[^ ■ e > |(^)* we conclude 

(^) * steps for the separation between M and 2M, provided \fM < 



that we need m < O 



[^yVN. The same technique will also work for the case K = aM , with some a > 1. 

Example 3 Note that if the conditions of the theorem hold for the triple M, K , N they 
also hold for uM , uK , uN where n is any integer. Also, the angles 6m, o-nd 6^ remain the 
same, and therefore so does the number of iterations required to complete the job, despite 
the fact that the database has increased n-fold. This means that we should take care only of 
triples M, K , N that do not have a common divisor. 

Example 4 Suppose that our information is that one of the following cases obtains: 
\S\ = Ml, or \S\ = M2, or ...,\S\ = M^. We can inductively use multiple Diophantine ap- 
proximations as in Eq ([^.- First find an odd integer I such that sm{l^^) ~ 1 for 1 < j < [|], 
while sin(/-^) ^ 0, for [|] < J < 1- If a measurement discovers a marked element then 
with high probability \S\ = Mj for some I < j < [|], otherwise it is one of the other cases. 
Now, divide the resulting set of possibilities into two halves and continue the process. After 
~ log2 r successive Diophantine approximations we are guaranteed to find the answer. The 
trouble is that the larger r is the larger I is likely to be, and it is not clear when the process 
yields better than classical outcomes. 

Example 5 If M = and K > 0, then we are just back with Grover's type algorithm. If 
\S\ = nothing happens to {ip) = \c()m' '^'^'^ ^/ 1*^1 = K we get close to 

Two remarks on the general problem should be made: Firstly, for arbitrary values of M 
and K the minimal size of the number of iterations m depends on our ability to obtain a 
lower bound on the uniform Diophantine approximation to the quotient 7 = By this we 
mean finding natural numbers p and q such that |p7 — q\ is small, but not too small. The 
reason will become clear from the proof below. Intuitively, if p^ stays close to an integer for 
a long segment of values of p, then the two numbers ^ and |^ become hard to separate with 
a small /. This means that the general separation problem runs into a difficulty similar to 
the problem of small denominators (or divisors) in classical mechanics . It is likely that 
a formulation of the algorithm in terms of a continuous adiabatic quantum computer will 



demonstrate more clearly the relation between our problem and the KAM-type of problems, 
in the sense that small divisors may show up in the spectral gap. Note also that some of 
these small divisor problems may be overcome by using the trick in Example 2, namely by 
adding artificial elements to the database and changing the values of Ok and 6m- 

Secondly, a remark about the actual value of / = 1{M, K, N), the number of iterations 
needed to complete the task. The proof below is giving a pretty good estimation of /. 
However, note that this is essentially a different problem. Once a "table" of the values 
of / is generated for the appropriate M, K, N, it can serve all search problems, no matter 
what the nature of the objects in the database, and the character of the oracle. Such 
"table" may allow us to decide what is the best strategy to use. We shall just briefly 
indicate how to formulate this problem algebraically: Denote Ti{x) = cos[/ cos~^(x)], then 
Ti is the / degree Chebyshev's polynomial (of the first kind) [11]. Using Eq Q we see that 
cos(/^) = Ti ( ^ ^^-j^) and similarly cos(/^) = Ti{^J^^-). To get rid of the square roots 
we can use the identity 2Tf{x) — 1 = T2i{x) = T/(2x^ ~ 1); so that finally our task is to find 
the smallest odd / such that Ti{^-^^^) ^ +1 while Ti{ ^~^^ ) ^ — 1. This observation may 
also assist in generalizing our result to other values of M , and A^. 



Proof of the theorem : Denote 7 = I77 = ""tt^^^P"? ^6 take the following three steps 



Step 1: If M < a: < f then 



< 7 - 1 < V2{^K/M - 1) (9) 

That 7 > 1 is obvious since M < K and sin~^(a;) is increasing. For the right hand 
estimation we use the mean value theorem. First note that if < 0:2 < 2:1 < 1 then 



sin ^(xi) ]^ _|_ ^(xi) — sin ^(x2) ^ ^ ^/T^^~x^ xi — X2 
sin"^(x2) sin"^(a;2) - sin~^(0) a/1 - xj ^2 

for some X3 and X4 such that Xi > X4 > X2 > X3 > 0. Now, substitute Xi = K/N and 

TV 
2 



X2 = and remember that 1 — Xg < 1, and 0:4 < ^, so that the condition K < 

entails 7 < 1 + V2{y/K/M - 1). 

Step 2: Choose p to be the nearest odd integer to ^^^^-^^ , then p — < 1, and also 

P < ^ (assuming 7 — 1 < 7). and altogether: 



Pl-P- ^ 



7 



Now, add the condition \l jj < 16(7 — 1)^. Denote by s the nearest odd integer to J^, 



then 



s — 



47r 
0M 



< 1, and put / = ps therefore I is also odd. Then, since sin ^[x) < for 
< a; < 1, we have from Eq pO|) 

1 



^-4 



< 



9k 


1 






Pit 









< 2(7-1) 



(11) 



47r 



P 



47r 



< (7 - 1) 



Step 3: Let £ > 0. To complete the proof all we have to do is impose the condition 
2(7 -I) <s. But then by Eq Q this will be satisfied if ^/K < (1 + ^)VM. To estimate 



/ note that by our definition p 



1 



4(7-1) 



while s 



4^ where 



indicates equality up to 



ps 



. Using once more the mean value theorem for sin ^(x) we get 



±1. Hence, I 

Returning to the issue of small denominators, consider how it is avoided in our proof: 
On the one hand 7 — 1 is small, indeed 7 — 1 < | is our basic constraint. On the other 
hand p is of the order of magnitude of (7 — 1)"^, and / = > p, so that 7 — 1 cannot be 
too small. This is the balance that should be struck if we wish to generalize the result to 
other values of 7 = we have to obtain a uniform Diophantine approximation |p7 — q\ 
which is small, but reasonably bounded from below. A further complication is that p has 
to be odd. General lower bounds of this kind exist for algebraic numbers, but 7 is typically 
transcendental. However, we are dealing with a very special case for which a good lower 
bound may exist. Also, we can move from a bad case to a better one by adding artificial 
elements to the database, as in Example 2. 

Conclusion Given an oracle that identifies the elements of a subset S C {1, 2, A^}, 
and knowledge that either \S\ = M ox \S\ = K ^ for M < K, we demonstrated how to 
decide which is the case by iterating Grover's rotation m < ^^"^^^ times. The algorithm is 
working for a certain range of values M, K, and A^, and employs simultaneous Diophantine 
approximations. This means that we conceive of Grover's algorithm as an orbit of an ergodic 
automorphism of the torus T^, and ask how quickly it enters a given open subset of T^. We 
showed how to apply this process in some special cases, and noted that in other cases the 
algorithm may be frustrated because of a 'small divisor' type of problem. 
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